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Abstract. We have computed new solar models using the 
same stellar evolution code as described in Charbonnel, 
Vauclair and Zahn (1992). This code, originating from 
Geneva, now includes the computation of element segrega- 
tion for helium and 12 heavier isotopes. It may also include 
any type of mixing of the stellar gas, provided this mix- 
ing can be parametrized with an effective diffusion coeffi- 
cient as a function of radius. Here we introduced rotation- 
induced mixing as prescribed by Zahn (1992). We present 
five solar models: 1) the standard model, computed with 
heavy element abundances as given by Grevesse (1991); 
2) a model including pure element segregation (no mix- 
ing outside the convective zone) with Grevesse (1991) 
as initial abundances; 3) same model as (2), but iter- 
ated so that the final abundances are those of Grevesse 
(1991); 4) a model with both element segregation and 
rotation-induced mixing, leading to lithium and beryllium 
depletion consistent with the observations, with Grevesse 
(1991) as initial abundances; 5) same model as (4) but 
iterated to obtain Grevesse (1991) as final abundances. 
This model (5) now represents our best new solar model 
consistent with the observations. 

The u = — function computed as a function of radius 
in these new solar models are compared to the helioseismo- 
logical results obtained for the same function by Dziem- 
bowski et al (1994). Improving the physics of the models 
leads to a better consistency with helioseismology. In our 
best model (5), which includes both segregation and mix- 
ing, the relative difference in the u function between the 
model and the helioseismological results is smaller than 
0.5 per cent at all radii except at the center and the sur- 
face. Meanwhile lithium is depleted by a factor 155 and 
beryllium by a factor 2.9, which is consistent with the 
observations. The bottom of the convective zone lies at 
a fractional radius of 0.716, consistent with helioseismol- 
ogy. The neutrino fluxes are not decreased in any of these 
models. 



The models including the computations of element seg- 
regation lead to a present surface helium abundance of: 
Ysurf between 0.248 and 0.258, which is in satisfactory 
agreement with the value derived from helioseismology. 

Key words: Physical data and processes : diffusion - Sun 
: abundances - evolution - interior - rotation 



1. Introduction 

The Sun is by far the most well known of all the stars. Its 
mass, radius, luminosity and age have been determined 
with a high degree of precision (Table 1). The mass is ob- 
tained from the motion of planets, the radius from eclipses 
(the value given in table 1 is reduced to an optical depth 
t=2/3), the luminosity from measurements of the solar 
constant above the earth's atmosphere. Some discussion 
remains about the solar age: it is generally taken as 4.6 
billion years although Guenther (1989) and Demarque and 
Guenther (1991) suggest a smaller age consistent with the 
oldest meteorites. The photospheric solar element abun- 



Table 1. Solar parameters 



mass 


(1.9891 ±0.00 04) x 10 33 g 


radius 


(6.959 ±0.001) x 10 10 cm 


luminosity 


(3.851 ± 0.005) x 10 33 ergs.s -1 


age 


4.6 ±0.15 Gyr 


mass loss 


2 x 10" 14 MQ.yr" 1 



dances are now precisely known, after the studies by An- 
ders and Grevesse (1989) modified by Grevesse (1991). See 
also Grevesse and Anders (1991). Lithium and beryllium 
are both depleted in the Sun with ratios : 
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with an relative uncertainty of 30% in both cases. Any 
consistent solar model must account for these depletion 
factors. 

The study of the internal structure of the Sun entered a 
new age with the birth of helioseismology. Millions of solar 
p-modes have been detected (including the (21 + 1) mul- 
tiplcts). An inversion of the measured frequencies yields 
accurate and detailed information about such structural 
functions as pressure, p (r), and density, p (r), in the Sun's 
interior. A particularly high precision is achieved in the 
determination of u = 2 throughout the Sun, and in the 
localization of the bottom of the convective envelope. No 
assumption regarding the transport of energy and chemi- 
cal elements is introduced at this stage of seismic sound- 
ing. The only essential assumption is that of mechanical 
equilibrium, which is partially testable by means of helio- 
seismology. 

Christensen-Dalsgaard et al. (1993) have shown that 
gravitational settling and element mixing processes signif- 
icantly affect the calculated speed of sound near the bot- 
tom of the convective envelope. These processes lead to a 
lower surface helium abundance - a quantity which may 
also be directly inferred form hclioscismic data if the equa- 
tion of state is specified. The results of many independent 
inversions clearly demonstrate that helium settling must 
take place in the Sun. A question we ask in this paper 
is whether helioseismology provides useful constraints on 
the mixing processes. 

An important result of helioseismology is the precise 
determination of the bottom of the solar convective zone: 
^ = 0.713 ± 0.003 (Christensen- Daalsgard et al. 
1993). This value corresponds to that obtained from the 
Schwarzschild criterium, leading to a strong constraint on 
overshooting (section 3). In the following we show that, 
contrary to a common idea generally spread among solar 
physicists, this constraint on overshooting is not a prob- 
lem for the explanation of the lithium depletion in the Sun. 
In any case, explaining the lithium depletion in the Sun 
by overshooting would not be consistent with the lithium 
observations in other stars (e.g. galactic clusters). 

Charbonnel, Vauclair and Zahn (1992) (CVZ) and 
Charbonnel et al. (1994) showed that the lithium defi- 
ciency in solar type stars can be accounted for by rotation- 
induced mixing. We will show here that such a mixing can 
also account for Li and Be depletion in the Sun, without 
destroying the consistency with helioseismology. 

The recipe to construct solar models is well known 
and has been explained many times in the literature (Cox, 
Guzik and Kidman 1989, Bahcall and Pinsonneault 1992, 
Turck-Chieze and Lopes 1993, Proffitt 1994 etc.). We be- 
gin with a homogeneous 1M model on the zero-age main 
sequence. Then the model is evolved step by step, by tak- 
ing into account the modification of the chemical com- 
position due to nuclear reactions. At the age of the Sun, 
it must reproduce the data given in table 1 within the 
observed uncertainties. To obtain such a high precision, 



two parameters are traditionally adjusted, with an itera- 
tion prodedure: the 4 He mass fraction Y and the ratio of 
the convective mixing length to the pressure scale height, 
generally referred to as a. 

The computation of the standard solar models includes 
the assumption that the stellar gas as a whole is in hy- 
drostatic equilibrium. This fundamental "first equation" 
of the internal structure of stars assumes that a blob of 
stellar gas is in equilibrium due to the effect of gravity 
downwards and pressure gradient upwards, which is cor- 
rect in first approximation. It does not take however into 
account the fact that the stellar gas is a mixture of many 
different species, which do not have the same weight. This 
is the problem of element segregation, which indeed is a 
fundamental process inherent to the stellar structure. 

In this paper we present new solar models computed 
with a stellar evolution code including element segregation 
and mixing, as described in CVZ. The physics included in 
these computations is discussed in section 2, a discussion 
about helioseismology appears in section 3 and the results 
are given in section 4. 

Five solar models will be discussed: 1) the best stan- 
dard model 2) a model including pure clement segregation 
(no mixing outside the convective zone) with initial abun- 
dances corresponding to the Grevesse (1991) mixture; 3) 
same model as (2), but iterated so that the final abun- 
dances are those of Grevesse (1991); 4) a model with both 
element segregation and rotation-induced mixing, leading 
to lithium and beryllium depiction consistent with the ob- 
servations, with Grevesse (1991) as initial abundances; 5) 
same model as (4) but iterated to obtain Grevesse (1991) 
as final abundances. This model (5) now represents our 
best new solar model consistent with the observations. 

All these models are compared to the helioseismolog- 
ical results obtained by Dziembowski et al (1994). It is 
very encouraging to see that improving the physics of the 
models leads to a very good consistency with helioseismol- 
ogy. In our best model (5), which includes both segrega- 
tion and mixing, the relative difference in the u function 
between the model and the helioseismological results is 
smaller than 0.5 per cent at all radii except at the cen- 
ter and the surface. Meanwhile lithium is depleted by a 
factor 150 and beryllium by a factor 2.9, which is consis- 
tent with the observations. The bottom of the convective 
zone lies at a fractional radius of 0.716, consistent with 
helioseismology. 

The computed neutrino fluxes will be presented in the 
tables of results for each models. They are not decreased 
in any of them. These results will be shortly discussed but 
a complete discussion of the solar neutrino problem is out 
of the scope of the present paper. 
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Table 2. Initial parameter and main physical surface parameter of solar models. 





Y 


a 


^surface 


X sur face 


L 

(lO^crg.s- 1 ) 


R 

(10 n cm) 


Li/Lio 


Be/Beo 


Model 1 


0.2782 


1.652 


0.2782 


0.7028 


0.385145 


0.695976 


1 


1 


Model 2 


0.2762 


1.776 


0.2477 


0.7341 


0.385154 


0.696368 


1/2.89 


1/1.17 


Model 3 


0.2798 


1.789 


0.2513 


0.7297 


0.385143 


0.695982 


1/3.50 


1/1.17 


Model 4 


0.2770 


1.761 


0.2563 


0.7252 


0.385131 


0.695980 


1/124.58 


1/2.88 


Model 5 


0.2793 


1.768 


0.2584 


0.7226 


0.384993 


0.695849 


1/155.03 


1/2.91 



Table 3. Main physical parameter of solar models at the base of the convective zone and at the center. 





r C z 
«0 




pcz 


Y c 


x c 


T c 


Pc 


Pc 




(10 6 K) 


(g.cm" 3 ) 






(10 6 K) 


(g.cm -3 ) 


( dyn.cm ~ 2 ) 


Model 1 


0.725 


2.100 


0.166 


0.6346 


0.3459 


15.56 


150.66 


2.303 10 17 


Model 2 


0.716 


2.158 


0.185 


0.6416 


0.3383 


15.63 


153.81 


2.344 10 17 


Model 3 


0.714 


2.178 


0.189 


0.6464 


0.3326 


15.70 


154.17 


2.345 10 17 


Model 4 


0.717 


2.162 


0.185 


0.6431 


0.3368 


15.63 


154.17 


2.350 10 17 


Model 5 


0.716 


2.175 


0.188 


0.6465 


0.3328 


15.67 


154.53 


2.350 10 17 



2. The computations 

2.1. Input microphysics 

- Equation of state : The present solar models are com- 
puted with the equation of state developped by Hummer 
& Mihalas (1988), Mihalas et al. (1988), Dappen et al. 
(1988), hereafter MHD. The MHD equation of state is 
based on the free-energy minimization method which im- 
plies that it is thermodynamically consistent. It treats 
pressure ionization carefully and takes into account non 
ideal effects such as Coulomb correction to pressure, pres- 
sure due to partially degenerate electrons and correction 
for size of particles. It also includes a large number of 
atomic, ionic and molecular species, with detailed parti- 
tion functions, containing weighted occupation probabili- 
ties. As shown by Christensen-Dalsgaard et al. (1988), the 
MHD equation of state highly improves the solar models 
in the helioseismological context. 

- Opacities : We use the OPAL radiative opacities by 
Iglesias et al. (1992) which include the spin-orbit inter- 
actions for Fe and relative metal abundances based on 
Grevesse (1991). These tables are complemented at low 
temperatures below 10 000K with the atomic and molec- 
ular opacities by Kurucz (1991). 

- Nuclear reactions : For hydrogen-burning we con- 
sider the tree pp chains and the CNO tri-cycle. We use the 
thermonuclear reaction rates Caughlan & Fowler (1988). 
Screening factors for the reaction rates are taken into ac- 
count according to the analytical prescription by Graboske 
et al. (1973). 



The key reaction for the generation of high energy neu- 
trinos and for the theoretical neutrino flux to be compared 
to the chlorine experiments results is 7 Be(p,7) 8 B. Caugh- 
lan & Fowler (1988) give a value of 0.0240 keV-barns for 
the corresponding low energy cross section factor, Si7(0), 
which is extrapolated from experimentally measured reac- 
tion cross sections. 

2.2. Element Segregation 

The process of element segregation in stars (also referred 
to as "microscopic diffusion") represents a basic physical 
process inherent to the stellar structure. As soon as the 
stars form out of gas clouds, they built density, pressure 
and temperature gradients throughout. Under such con- 
ditions, the various chemical species present in the stellar 
gas move with respect to one another, unless macroscopic 
motions force the chemical homogeneization. 

Although recognized by the pioneers of the study of 
stellar structure (Eddington, 1916 and 1926, Chapman, 
1917), this process was long forgotten in the computa- 
tions of stellar models, except for white dwarfs (Schatz- 
man, 1945). Only with the discovery of large abundance 
anomalies in main-sequence type stars (the so-called Ap 
and Am stars), which present characteristic variations 
of chemical elements with the effective temperature, was 
microscopic diffusion brought into light fifty years later 
(Michaud, 1970, see other references in Vauclair and Vau- 
clair, 1982). 

At that time, the effects of microscopic diffusion were 
supposed to be important only when the diffusion time 
scale was smaller than the stellar age. In the Vauclair and 
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Fig. 1. Standard Solar model (model 1). The top graph rep- 
resents the difference between the u function (u — —) deduced 
from helioseismology and the computed one. The ordinates rep- 
u (seismic) -« (model) 
u (seismic) 

The bottom graph shows the abundance profiles for some ele- 
ments in the two models (The mass fractions are multiplied by 
250 for 3 He, by 200 for 12 C, by 1000 for 13 C, by 140 for 14 N 
and by 75 for 16 O). 



Vauclair (1982) review paper, Fig. 1 shows the regions in 
the HR diagram where microscopic diffusion could lead 
to "observable" abundance variations. The Sun was ex- 
cluded, although at the border of the "permitted domain" . 
In the present days, due to helioseismology, abundance 
variations of the order of a few percent become indirectly 
detectable: we have entered a new area in this respect. 

Several authors have computed the gravitational and 
thermal diffusion of helium and heavier elements in 
the Sun with various approximations (see references in 
Michaud and Vauclair, 1991). More recently the influence 
of diffusion on the solar oscillation modes and on the solar 
neutrino fluxes have been studied in various ways (Cox, 
Guzik, Kidman, 1989; Bahcall and Pinsonneault, 1992; 
Proffitt, 1994; Thoul, Bahcall and Loeb, 1994). 

Element segregation represents in fact a competition 
between two kinds of processes. First the atoms move un- 
der the influence of external forces (due to gravity, ra- 
diation, etc.), second they collide with other atoms and 
share the acquired momentum with them in a random way, 
which slows down their motion. This competition leads to 
element stratification and decreases the entropy. 

The computations of microscopic diffusion are based 
on the Boltzmann equation for dilute collision-dominated 
plasmas. At equilibrium the solution of the equation is the 
maxwellian distribution function: / = /(0). In stars the 
distribution is not maxwellian, but the deviations from 
the maxwellian distribution are very small. Two differ- 
ent methods arc devclopped to solve the Boltzmann equa- 
tion in this framework. The first method relies on the 
Chapman-Enskog procedure (described in Chapman and 
Cowling, 1970), with convergent series of / computed with 
successive approximations. This method is used, for ex- 
ample, by Bahcall and Loeb (1990), Proffitt and Michaud 
(1991), Michaud and Vauclair (1991), Bahcall and Pinson- 
neault (1992) and CVZ, in which a complete description 
of the numerical schemes may be found. 

The second method is that of Burgers (1969), in which 
separate flow and heat equations for each component of a 
multi-component mixture are solved simultaneously. De- 
scriptions of this method may be found for example in 
Cox, Guzik and Kidman (1989), Proffitt and Vanden- 
Berg (1991), Thoul, Bahcall and Loeb (1994), Richer 
and Michaud (1993). This method does not include for 
the moment the problem of partial ionisation which has 
been studied within the framework of the Chapman and 
Cowling method (Montmerle and Michaud 1976; Vauclair, 
Hardorp, Peterson, 1979; Alecian and Vauclair, 1981). 

In the present paper the first method has been used for 
the treatment of the Boltzmann equation, as described in 
CVZ. For collisions between charged ions, the Paquette et 
al. (1986) method has been introduced. The basic question 
concerns the divergence of the coulomb interaction cross 
sections. In the first computations of diffusion, the "Chap- 
man and Cowling approximation" was used, assuming a 
cut-off of the cross section equal to the Debye shielding 
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Fig. 2. Models with pure element segregation. The presentation is the same as in figure 1. The left graphs represent the model 
obtained with initial abundances as given by Grevesse (1991) (model 2). The right graphs represent an interated model computed 
so that the final abundances are very close to Grevesse's observed ones (table 4). The Ati s are clearly smaller for the iterated 
model. The abundance variations are shown in the bottom graphs. The relative abundance variations in the convective zone 
compared to the interior (depletion of all the elements except hydrogen) are the same in figures 2b and 2d, but the absolute 
values differ, so that in fig. 2d the outer values are close to those given in fig. lb for all elements except helium 4, which is 
adjusted in the calibration of the model. 
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Table 4. Initial abondances, final surface abondances, and Grevesse mixture in mass fraction. 



Model 2 Model 3 





Initial 


Grevesse 


Final 


Grevesse 


Initial 


Grevesse 


Final 


Grevesse 




mixture 


(1991) 


mixture 


(1991) 


mixture 


(1991) 


mixture 


(1991) 




used 


with same 


obtained 


with same 


used 


with same 


obtained 


with same 




(12) 


Y than 12 


(F2) 


Y than F2 


(13) 


Y than 13 


(F3) 


Y than F3 


X 


0.704800 


0.704800 


0.734103 


0.733309 


0.700428 


0.701185 


0.729746 


0.729653 




0.000042 


0.000042 


0.000040 


0.000044 


0.000046 


0.000042 


0.000043 


0.000044 


Y 


0.276200 


0.276200 


0.247691 


0.247691 


0.279815 


0.279815 


0.251347 


0.251347 


X C 12 


0.003288 


0.003288 


0.003084 


0.003421 


0.003593 


0.003271 


0.003373 


0.003404 


X C 13 


0.000057 


0.000057 


0.000053 


0.000059 


0.000063 


0.000057 


0.000059 


0.000059 


Xj\Tl4 


0.000976 


0.000976 


0.000919 


0.001015 


0.001062 


0.000971 


0.001001 


0.001010 


Xjyis 


0.000004 


0.000004 


0.000004 


0.000004 


0.000004 


0.000004 


0.000004 


0.000004 


X 16 


0.009494 


0.009494 


0.008971 


0.009878 


0.010313 


0.009445 


0.009754 


0.009828 


X G 17 


0.000007 


0.000007 


0.000007 


0.000007 


0.000007 


0.000007 


0.000007 


0.000007 


XqIS 


0.000023 


0.000023 


0.000022 


0.000024 


0.000025 


0.000023 


0.000023 


0.000024 


X Afe 20 


0.001603 


0.001603 


0.001603 


0.001668 


0.001657 


0.001595 


0.001657 


0.001660 


X Afe 22 


0.000140 


0.000140 


0.000140 


0.000146 


0.000145 


0.000140 


0.000145 


0.000145 


X Mg 24 


0.000497 


0.000497 


0.000497 


0.000517 


0.000514 


0.000494 


0.000514 


0.000515 


X Mg 25 


0.000068 


0.000068 


0.000068 


0.000071 


0.000071 


0.000068 


0.000071 


0.000071 


X M g26 


0.000081 


0.000081 


0.000081 


0.000085 


0.000084 


0.000081 


0.000084 


0.000084 


Z/X 


0.026958 


0.026958 


0.024800 


0.025910 


0.028207 


0.027097 


0.025909 


0.026040 



Model 4 Model 5 





Initial 


Grevesse 


Final 


Grevesse 


Initial 


Grevesse 


Final 


Grevesse 




mixture 


(1991) 


mixture 


(1991) 


mixture 


(1991) 


mixture 


(1991) 




used 


with same 


obtained 


with same 


used 


with same 


obtained 


with same 




(12) 


Y than 12 


(F2) 


Y than F2 


(13) 


Y than 13 


(F3) 


Y than F3 


X 


0.704016 


0.704016 


0.725216 


0.724734 


0.701229 


0.701714 


0.722576 


0.722551 


X ffe 3 


0.000042 


0.000042 


0.000076 


0.000043 


0.000018 


0.000042 


0.000054 


0.000043 


Y 


0.276984 


0.276984 


0.256266 


0.256266 


0.279286 


0.279286 


0.258449 


0.258449 


X C 12 


0.003284 


0.003284 


0.003149 


0.003381 


0.003498 


0.003273 


0.003355 


0.003371 


X C 13 


0.000057 


0.000057 


0.000054 


0.000059 


0.000062 


0.000057 


0.000059 


0.000059 


Xjvl4 


0.000975 


0.000975 


0.000938 


0.001003 


0.001035 


0.000971 


0.000996 


0.001000 


Xjyis 


0.000004 


0.000004 


0.000004 


0.000004 


0.000004 


0.000004 


0.000004 


0.000004 


X D 16 


0.009483 


0.009483 


0.009142 


0.009762 


0.010052 


0.009452 


0.009693 


0.009733 


X Q 17 


0.000007 


0.000007 


0.000007 


0.000007 


0.000007 


0.000007 


0.000007 


0.000007 


XqIS 


0.000023 


0.000023 


0.000022 


0.000024 


0.000025 


0.000023 


0.000024 


0.000024 


X Afe 20 


0.001602 


0.001602 


0.001602 


0.001649 


0.001641 


0.001596 


0.001641 


0.001644 


X Afe 22 


0.000140 


0.000140 


0.000140 


0.000144 


0.000144 


0.000140 


0.000144 


0.000144 


X Mg 24 


0.000496 


0.000496 


0.000496 


0.000511 


0.000510 


0.000495 


0.000510 


0.000510 


X Mg 25 


0.000068 


0.000068 


0.000068 


0.000070 


0.000070 


0.000068 


0.000070 


0.000070 


X Mg 26 


0.000081 


0.000081 


0.000081 


0.000084 


0.000084 


0.000081 


0.000084 


0.000083 


Z/X 


0.026988 


0.026988 


0.025534 


0.026217 


0.027787 


0.027077 


0.026260 


0.026296 



length. Paquette et al. (1986) proposed a more precise 
treatment of this problem, with a screened coulomb po- 
tential in which the characteristic length is taken as the 
largest of the Debye length and interionic distance. The 
Paquette et al. tables of collision integrals have been ex- 
tensively used in the present paper. 

The radiative acceleration on the elements are not 
included in these computations. From crude approxima- 



tions, we suppose that they are negligible in solar type 
stars (Michaud et al., 1976), but this should be tested 
in the future, as more precise computations may lead to 
larger values (Michaud, 1987). 

2.3. Rotation-induced mixing 

We introduced in our computations the rotation-induced 
mixing as prescribed by Zahn (1992). In a rotating star, 



O. Richard et al.: New Solar Models Including Helioseismological Constraints and Light-Element Depletion 7 



Table 5. Predicted neutrino flux in 10 10 cm 2 .s 





<E>(pp) 


<E>(pep) 


$( 7 Be) 


$( 8 B) 

V / 




$( 13 N) 


$( 15 0) 


$( 17 F) 


Model 1 


6.00 


1.39 10~ 2 


0.45 


5.45 1(T 


-4 


4.52 10~ 2 


3.85 10" 2 


4.90 10 -4 


Model 2 


5.96 


1.39 10~ 2 


0.47 


5.98 10" 


■4 


5.05 10~ 2 


4.31 10~ 2 


5.53 10~ 4 


Model 3 


5.94 


1.38 10~ 2 


0.48 


6.38 10" 


■4 


5.79 10~ 2 


4.98 10~ 2 


6.40 10~ 4 


Model 4 


5.96 


1.39 10~ 2 


0.47 


6.06 10" 


■4 


5.09 10~ 2 


4.35 10~ 2 


5.60 10" 4 


Model 5 


5.94 


1.38 10~ 2 


0.48 


6.33 10" 


■4 


5.59 10~ 2 


4.81 10~ 2 


6.18 10~ 4 



Table 6. Predicted neutrino capture rates for the chlorine and gallium experiments in SNU. 



Neutrino Model 1 Model 2 Model 3 Model 4 Model 5 



source ($a)ci ($Q-) G a (ggjcj (^) Ga ($a)ci (^) Ga (gg)gi (®v) Ga ($a)ci (fropGa 



PP 





70.722 





70.353 





70.042 





70.303 





70.068 


pep 


0.223 


2.992 


0.223 


2.992 


0.221 


2.965 


0.223 


2.993 


0.221 


2.973 


7 Be 


1.064 


32.722 


1.117 


34.342 


1.149 


35.352 


1.124 


34.574 


1.146 


35.255 


8 B 


5.781 


13.252 


6.34 


14.534 


6.768 


15.515 


6.42 


14.718 


6.706 


15.372 


13 N 


0.075 


2.792 


0.084 


3.118 


0.096 


3.576 


0.084 


3.143 


0.093 


3.456 


IBq 


0.254 


4.46 


0.284 


4.994 


0.329 


5.775 


0.287 


5.045 


0.318 


5.581 




0.003 


0.057 


0.004 


0.065 


0.004 


0.075 


0.004 


0.066 


0.004 


0.072 


Total 


7.4 


126.997 


8.052 


130.398 


8.567 


133.3 


8.142 


130.842 


8.488 


132.769 



due to centrifugal effects, the gravity equipotentials are 
no more spherical, which induces a circulation of matter 
between polar and equatorial regions: the so-called merid- 
ional circulation. This circulation itself induces a transport 
of angular momentum, thereby creating shears which be- 
come unstable in the horizontal direction, while the ver- 
tical shears are stabilized by the density gradient. This 
large scale horizontal turbulence decays into small scales 
and becomes 3D when the turnover rate of the turbulence 
exceeds the angular velocity. 

Meanwhile the horizontal turbulence "cuts down" the 
effect of advection on the transport of the chemical species, 
as the elements which go up in the upward flow of matter 
can be transported into the downward flow by horizontal 
motions before reaching the top layers. The transport of 
angular momentum is more efficient than the transport of 
chemicals. In the limit of extremely large horizontal diffu- 
sivity, the chemical composition is constant along a level 
surface, and the transport of chemicals is negligible. The 
angular momentum behaves differently as, when the rota- 
tion velocity is constant along a level surface, the angular 
momentum is not, so that it is still transported. 

The horizontal transport of angular momentum 
smoothes out the original meridional circulation. Taking 
this feedback effect into account, Zahn (1992) showed that 
the whole process is stopped within an Eddington-Sweet 
time-scale, unless angular momentum is extracted from 
the star due to a wind. 

In case of a moderate wind which extracts angular mo- 
mentum at the rate (^) , an asymptotic regime is reached 
with a circulation velocity as a function of radius given by 



(following Zahn (1992) eq.4.15): 

5 



U(r) = 



?, fdJ\ 



p(r) r 4 tt 8tt \dt J 
where p is the local density, or: 

Mr) = - - ^ ( -) 

y ' 2 fi(r) r M(r) p(r) \dt J 



(1) 



(2) 



where f2(r) is the local angular velocity, M(r) the mass 
inside radius r and p m the average density inside radius 
r. 

The effective diffusion coefficient is then expressed by 
Zahn's cq. 4.21 : 



Deft = - 



C h r\U(r)\ 



50 



(3) 



where Ch is a parameter related to the horizontal viscosity 
{Ch Si 1) and a is related to the differential rotation: 



1 dln(r 2 n) 
a = — 

2 din r 



(4) 



If the deviation of from solid rotation is neglected 
(a = 1), and if is supposed to decrease with time fol- 
lowing a "Skumanich law" for which £1 oc we obtain: 



dJ 
dt 



cx 



and the effective mixing coefficient is of the form: 
Deff - r.U(r) cx — 



(5) 



(6) 
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Equation (3) has been used in the present computations 
with a = 1. The proportionality factor Cn has been ad- 
justed to obtain the right lithium depletion in the Sun 
(section 4). 

2-4- Stabilizing /i- gradients 

Mixing processes in stars may be stabilized in the re- 
gions where the mean molecular weight rapidly decreases 
with increasing radius. This occurs specially in the nuclear 
bruning core: we can thus infer that the rotation-induced 
mixing becomes inefficient as soon as the /^-gradient be- 
comes larger than some critical value. This question has 
been discussed by several authors (Mestel (1965), Hup- 
pert and Spiegel (1977)). Although no precise value can 
be given for this critical /i-gradient, an order of magnitude 
can be obtained from simple considerations. 

Huppert and Spiegel (1977) suggest that mixing can 
penetrate the nuclear burning core within a scale height 
given by: 

fi(r) 



h ~ r 



(7) 



where r is the local radius, £l(r) the angular rotation veloc- 
ity and the buoyancy frequency due to the ^-gradient. 



N 2 ~ 



GM{r) 



din n 



dr 



(8) 



From eq. (7) we can derive a critical /i-gradient ob- 
tained by specifying that h must be a small fraction of r 

(h = e r): 



(vw )c 4^ 2(rJ 



(9) 



e 2 GM(r c ) 

where all the quantities should be computed at the place 
where the actual V In [i is equal to the critical one. 

With values of r c between 0.1 Rq and 0.2 Rq and val- 
ues of n between 3 x 10~ 6 (V Rq ~ 2 km.s" 1 ) and 10~ 4 
(Vr ~ 70 km.s" 1 ) we find, for s = 0.1: 



4 x 10~ 15 < (Vln n) c < 4 x 10 



-12 



Solar structure computations lead to stronger constraints 
on the critical /z-gradient. With (Vln fj,) c <^ 10~ 14 no mix- 
ing could occur in the Sun after 0.1 Gyr. With (V In fi) c <; 
10~ 12 too much mixing would occur in the core, and 
the consistency with helioseismology would be lost (Gaige 
1994). Our best model is obtained with (Vln /x) c = 4 x 
10~ 13 (section 4.3). 



2.5. The diffusion routine 

The stellar evolution code used in these computations is 
the Geneva code, described several times in the literature. 
The system of nuclear reactions and the abundance vari- 
ations in the standard models are computed as in Maeder 
(1983). They are separatly determined for 15 isotopes : H, 
3 He, 4 He, 12 C, 13 C, 14 N, 15 N, 16 0, 17 0, 18 0, 20 Ne, 22 Ne, 



24 Mg, 25 Mg, 26 Mg. The heavier elements are combined in 
a single mass fraction Z. 

We have added in this code a diffusion routine for each 
isotope in a similar way as described in CVZ. The diffusion 
equations are written in lagragian coordinates as : 



dcj _ d 2 c t 
dt 1 dm 2 



dcj 
dm 2 



Fi 



(10) 



where a stands for the concentration of isotope i, and: 



A 
F 



(47T7-V) 2 (D eff + Dl) 

dr 2 p(D cff D* s ) _ 
dm 

^ dm 



(4irr 2 py [An 



Vi (11) 



Here D e f f is the effective mixing coefficient (the same 
for all the isotopes) while D\ is the segregation coefficient, 
computed for each isotope using Paquette et al. (1986) ta- 
bles. Vi represents the segregation (microscopic) velocity. 
The nuclear destruction rate Xi is only included in Fi for 
lithium and beryllium, which arc treated separately from 
the network. 

For the 15 isotopes included in the nuclear network, 
the computation procedure is the following: 

- at each evolutionary step, equation (10) is solved sep- 
arately for all the isotopes except hydrogen. The method is 
the same as described in CVZ: a Cranck-Nicholson scheme 
with the inversion of a tridiagonal matrix including all the 
mesh points down to the center 

- this diffusion routine is used with a smaller time step 
as the evolution time step, for a better precision. Typically 
20 resolutions of the diffusion equation are done between 
two computations of a complete model 

- the new mass fractions of each isotope are computed, 
taking into account the normalisation equation: 



X 



k^i 



Satisfying this equation needs a consistent resolution of 
the abundances of all the considered isotopes with the 
inversion of a (k + 1) order matrix 

- the abundance variations due to the nuclear reactions 
are then computed. 

3. Helioseismological constraints 

We shall compare the u (r) = P p £j function and the Y sur f 
quantity in our theoretical models with those in the seis- 
mic model derived from the observed p-mode frequencies 
with the method described by Dzicmbowski et al. (1994). 
The seismic model that we adopt here differs from the one 
presented in that paper in two respects. First we use new 
frequency data for I < 3 degrees from BISON network 
(Elsworth et al., 1994). This has a most noticeable effect 
in u (r) in the core. The second is the use of OPAL rather 
than MHD equation of state (EOS) in the reference model 
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Fig. 3. Models with element segregation and rotation-induced mixing. The presentation is the same as in figures 1 and 2. The 
left graphs correspond to models computed with initial Grevesse (1991) abundances (model 4). The right graphs correspond 
to iterated models so that the final abundances are close to the observed Grevesse (1991) abundances (model 5). Here the 
parametrisation of the rotation-induced mixing is that proposed by Zahn (1992) in the asymptotic regime. The mixing efficiency 
is supposed annihilated when the /i-gradient is larger than 4 x 10 -13 cm -1 ( see text). This cut-off is very important as a deeper 
mixing inside the nuclear burning core destroys the consistency of the computed and observed u values. Here the Am values are 
slightly smaller in figures 3a and 3c than in figures 2a and 2c. Meanwhile lithium and beryllium are destroyed as observed in 
the Sun (see figure 4). 
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adopted for the inversion (for consistency with the opaci- 
ties), which does not change the u values in a significant 
way. The only important consequence of this change is a 
somewhat higher seismic value for Y sur f. 

Within the adiabatic part of the convective zone the 
sound speed is determined solely by the value of the M to 
R ratio and by the value of adiabatic exponent, T\. Com- 
parison of the model and the seismic u (r) in this region 
provides therefore a test of the EOS. The situation in the 
radiative interior is more complicated. Let us note that 
approximately u oc — , and therefore the radiative trans- 
port of energy and the element diffusion directly affect 
u (r) below the convective envelope. 

In general, it is not possible to disentangle the un- 
certainty in the opacity and in the diffusion coefficients. 
Only the convective overshooting leading to a discontinu- 
ity in 4^ leaves a signature in the sound speed which is 
seismically detectable. Monteiro et al. (1994) looked for it 
and came up with an upper limit on the extent of such 
overshooting, which is 0.07 of the local pressure distance 
scale or 0.006 R . In our models such overshooting has 
been ignored and the above results show that this is a 
good approximation. The bottom of the convective zone in 
such case may be determined quite accurately. In the seis- 
mic model used in this work r cz = (0.7137 ± 0.0002) R Q , 
which agrees very well with the first accurate helioseismic 
value r cz = (0.713 ± 0.003) R Q (Christensen-Dalsgaard et 
al.,1993). 

There are still some uncertainties due to the opaci- 
ties. Unfortunately, we do not have a good way to assess 
limits for the induced modifications. Alternative opacities 
(OP, Seaton et al. 1994) are by about 40 percent lower at 
temperature and density ranges of interest here. However, 
this large difference may perhaps be due to a neglect of 
the plasma effects on the atomic properties, in the OP 
calculations, and therefore may not serve as an estimate 
of the uncertainty. 

In such a situation it is important to make use of the 
additional constraint which is the value of Y sur f deter- 
mined by means of helioseismology. In this comparison 
the uncertainty in the opacity data is rather unimportant. 
The possibility of a seismic measure of helium abundance 
rests on a large value of the Ti derivative with respect 
to Y in the Hell ionization zone and obviously requires 
very accurate thermodynamical data. The seismic values 
obtained with Dziembowski et al. (1994) method of inver- 
sion for two versions of the EOS adopted in the reference 
model are 

Y surf = 0.2440 ± 0.0003 for MHD 

Y surf = 0.2505 ± 0.0003 for OPAL 

The test of Ti in the lower convective zone points 
to OPAL data as more accurate, which should be ex- 
pected as the OPAL EOS is obtained in a more funda- 
mental way. We stress that the model values of Y sur f is 
quite insensitive to the choice of the EOS. The errors 
given above reflect only the frequency errors quoted by 
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Fig. 4. Lithium (top) and beryllium (bottom) variations pro- 
files in the best model 5 as a function of age. At the age of the 
Sun lithium is destroyed by 155 and beryllium by 2.9. 
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the observers. The actual uncertainty in the seismic val- 
ues is much larger. For instance, Basu and Antia (1995) 
find 0.246 and 0.249 for the corresponding quantities. The 
results are indeed depended on the adopted method of 
inversion. In particular they depend on smoothing the 

— function which describes the relative differences be- 
lt 

tween the solar and the model u (r) . Without smoothing, 
which is perhaps a better choice is the only goal is the 
determination of the He abundance, our method yields 
Y surf = 0.2548 for the OPAL EOS. (Pamyatnykh, private 
communication). The problem certainly requires further 
examination because, as we shall see, it is essential to re- 
duce the uncertainty to the 10 -3 level. 

4. The results 

Five solar models have been computed and compared to 
the helioseismological sun (Figures 1 to 3). The values of 
the characteristic parameters of the models, abundances 
and neutrino production are given in tables 2 to 6. 

4-1. The "standard" model 

This model includes the physics as discussed in section 2- 
1, with no element segregation and no other mixing than 
inside the convection zone. The comparison with helioseis- 
mology is presented in figure la, and the fractional abun- 
dance of several interesting elements in figure lb. The two 
curves in figure la correspond to the uncertainty in the 
helioseismological inversion. 

A comparison of seismic values of r cz and Y sur f with 
those for Model 1, given in Tables 2 and 3 reveals large 
differences. Also large differences are seen in the sound 
speed behavior throughout the whole interior as seen in 
Fig. 1. We do not pay attention to the differences in the 
outer part of the convective zone. They may be accounted 
for by inadequacies of the mixing length theory (outermost 
part) and in the MHD equation of state. Our main concern 
is the value ^ of about 0.015 in the outer part of the 
radiative interior, which is too large. 

Of course lithium is not depleted in model 1, neither 
is beryllium, and the solar neutrino fluxes are too high 
compared with the results of the solar neutrino detectors 
(tables 5 and 6). 

4-2. Models with pure element segregation (no mixing) 

- the first of these models (model 2) has been obtained 
with element segregation of helium and the initial abun- 
dances as given by Grevesse (1991). 

We see at once, from figure 2a, that the discrepancy 
below the convective zone is considerably reduced com- 
pared to model 1, while figure 2b shows the influence of 
segregation on the fractional abundances of the elements. 
The basic reason for the improvement in the u value is due 
to helium diffusion, which leads to a diminution of the he- 
lium abundance of about 10% inside the convection zone 



and below, thereby decreasing the local mean molecular 
weight. 

- model 2 is however inconsistent with the observed 
present element abundances due to the segregation. Model 
3 is similar to model 2, but iterated so that the final abun- 
dances correspond to the values given by Grevesse (1991). 
Table 4 shows the final abundances obtained, compared to 
the observations. It is also interesting to compare figure 
2d with figure 2b to see the difference in the fractional 
abundances of the elements. 

There is a dramatic improvement in the agreement 
with the seismic model once the element segregation is 
introduced. In Model 3 the coincidence of r cz and Y sur f 
with the seismic values is so perfect that, admittedly, it 
may be to some extent coincidental. 

From tables 2 to 6 we check that models 2 and 3 
are well calibrated, while the neutrino fluxes are not de- 
creased, as expected. 

We insist on the fact that these models are computed 
without any arbitrary parameters: element segregation is 
a simple consequence of the stellar physics, with no special 
assumption added. They should not be considered as non- 
standard models, but as improved standard solar models. 

Lithium and beryllium are not more depleted than he- 
lium in these models, as the depiction is only due to cle- 
ment segregation. It is necessary to introduce a mild mix- 
ing below the convection zone to account for the observed 
abundances of these elements. 

4-3. Models with element segregation and rotation-induced 
mixing 

Rotation-induced mixing as described in section 2 have 
been added in the computations of models 4 and 5. Here 
a parameter has been ajusted (namely the Ch factor in 
Zahn's prescription) so as to obtain a lithium depiction 
as observed. In consequence these two models must be 
considered as non-standard. 

In both models the mixing coefficient has been cut off 
at the core for a ^-gradient of 4 x 10~ 13 . Models with 
mixing down to the solar core lead to large discrepancies 
with helioseismology. On the other hand, in our models, 
the mixing is inefficient in the region of energy production 
in the Sun, and the consistency with helioseismology is 
even better than for models without mixing (see figures 3a 
and 3c). 

Model 4 is obtained with initial abundances as in 
Grevesse (1991), while model 5 has been iterated so that 
the final abundances are those of Grevesse (1991) (see ta- 
ble 4). It is interesting to note that model 5 is still closer 
to the helioseismological Sun that model 4: in all cases 
improving the physics leads to better results compared to 
the observations, which is very encouraging. 

While an inclusion of the element mixing somewhat 
improves the agreement in the u (r), it slightly destroys 
the perfect coincidence in r cz and Y sur f. We do not con- 



12 



O. Richard et al.: New Solar Models Including Helioseismological Constraints and Light-Element Depletion 




Fig. 5. Comparison between the parameters log T, log P, /i and log p in our "best model" 5 (solid lines) and the standard 
model (dashed lines). While no clear difference can be seen at these scales for T, P and p, the difference in jj, due to the element 
segregation is clearly visible. It is the basic reason for the better agreement of this model with the helioseismological sound 
speed, compared to the standard model. 
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sider any of these differences as really significant. Let us 
note that the seismic value of Y sur f obtained without reg- 
ularization, which we regard as more credible, is right in 
the middle of the values for model 3 and 5. The difference 
between the two models is 0.005. An accuracy of 0.001 in 
seismic Y sur f would thus yield an interesting constraint 
on the element mixing. The data certainly allow such ac- 
curacy. The problems lie in data analysis and in reliability 
of the EOS. 

In model 5 lithium is depleted with a ratio 1/155 and 
beryllium with a ratio 1/2.9, which is very close to the ob- 
servations. The lithium and beryllium abundance profiles 
for various times are given in figure 4. 

The difference in /x- values from models 1 and 5 is given 
in figure 5. 

5. Conclusion 

Model 5 is presently our best model. It is a well calibrated 
model, including element segregation and a parametrized 
rotation-induced mixing. It leads to a very good fit be- 
tween the computed u = ^ function and that deduced 
from hclioscismology below the solar convection zone. The 
consistency is also very good inside the convection zone ex- 
cept at the surface where a better mixing treatment should 
be introduced. Lithium and beryllium are depleted as ob- 
served. The neutrino fluxes remain too high. 

Below a radius of 0.4 of the solar radius a discrepancy 
remains in the u curves. Although the uncertainty on u is 
much larger in the core than at the surface, these features 
seem significant. This specific problem will be addressed 
in a forthcoming paper. 

Element segregation certainly takes place in the Sun's 
interior. It involves no nonstandard physics and its oc- 
curence is fully confirmed by the results of helioseismic 
inversion. On the other hand the macroscopic mixing is a 
hypothetical effect which provides a natural explanation 
of the Li and Be deficit but demands convincing observa- 
tional confirmations and constraints. At present stage we 
may only conclude that there is no conflict between mod- 
els reproducing solar Li and Be abundances and helioseis- 
mic data. There are good prospects for obtaining stringent 
constraints on the hypothetical mixing processes from hc- 
lioscismology. On the road to this goal we regard as most 
important to improve reliability and accuracy of the He 
abundance determination. Progress should also be made 
in assessing uncertainties in the opacity data, so that the 
information about the sound speed behavior may lead to 
a more direct probe of the He distribution in the outer 
part of the radiative interior. 

Note: Model 5 is available on request by electronic mail 
at the address: richard@obs-mip.fr. 
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